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Abstract 

We study the geometric phase acquired by an inertial atom whose trajectories are parallel to a 
reflecting boundary due its coupling to vacuum fluctuations of electromagnetic flelds, by treating 
the atom as an open quantum system in a bath of the fluctuating vacuum fields, and show that 
the phase is position dependent as a result of the presence of the boundary which modifies the 
field quantum fiuctuations. Our result therefore suggests a possible way of detecting vacuum 
fiuctuations in experiments involving geometric phase. 

PACS numbers: 
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Quantum theory has profoundly changed our conception of vacuum as a synonym of 
nothingness. As an inevitable consequence necessitated by the uncertainty principle, vac- 
uum fluctuates and thus may have rich structures. An intriguing issue is what are the 
physical consequences of vacuum fluctuations that exist all the time and whether these fluc- 
tuations can be directly detected. In this regard, let us note that the effects of vacuum 
fluctuations in free space may not be always observable, let alone a direct observation, since 
some physical quantities, energy for instance, are not well-defined in vacuum and one has to 
invoke certain renormalization schemes to make them finite. However, changes in the vac- 
uum fluctuations, e.g., those caused by the presence of boundaries, usually exhibit normal 
behaviors and can produce observable effects. The Lamb shift [1] and the Casimir [2] (and 
the Casimir- Polder [3, 4]) effects are the two most well-known examples, which have been 
precisely measured. Other examples of novel effects that arise as a result of the vacuum 
fluctuations include but are by no means limited to the light-cone fluctuations when gravity 
is quantized [5], the Brownian motion of test particles in an electromagnetic vacuum [6], 
and modifications of radiative properties of atoms in cavities such as the natural lifetimes 
and energy level shifts which have been demonstrated in experiments [7]. In this paper, we 
show that the vacuum fluctuations may also be directly detected though the measurement 
of geometric phase. As a related issue of vacuum fluctuation detection, it is worth noting 
that the quantum vacuum fluctuation in the position of a mechanical system has recently 
been clearly detected using a nanomechanical resonator [8]. 

The geometric phase is an important concept in quantum theory. In 1984, Berry studied 
the dynamics of a closed quantum system whose Hamiltonian varies adiabatically in a cyclic 
way, and found that there is, besides the familiar dynamical phase, an additional phase 
due to the geometry of the path enclosed during the unitary evolution of the system in 
the parameter space [9]. Ever since the inception, the geometric phase has aroused broad 
interest and has been extensively studied, both theoretically and experimentally [10]. Recent 
concerns about the geometric phase mainly focus on its potential of performing fault-tolerant 
quantum computation [11]. Because of the inevitable interactions between the qubits and 
the environment, a pure state will generically be driven to a mixed state. To deal with the 
effect of the environment, many attempts have been made which generalize the geometric 
phase of a closed system undergoing a unitary evolution to an open quantum system which 
undergoes a nonunitary one [12-16]. Remarkably, experiments have demonstrated both the 
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geometric phase of a mixed state undergoing a cyclic unitary evolution [22, 23] and that of 
an open system undergoing nonunitary one [24]. 

In fact, the impact of environment on the geometric phase is a crucial issue in any prac- 
tical implementations of quantum computing. The effects of different kinds of decoherence 
sources on the geometric phase have been analyzed [17-21]. It is remarkable that when 
the temperature is absolute zero, there is still a correction to the geometric phase caused 
by the environment, i.e., a reservoir of vacuum fluctuations [19-21]. However, this kind of 
inevitable vacuum fluctuation induced geometric phase is in general unobservable, since any 
phase variation is observed usually via some kind of interferometry between the involved 
state and certain selected reference states which are both inseparably and equally coupled 
to vacuum. Nevertheless, if, somehow, vacuum fluctuations are modified, then the geometric 
phase of the nonunitary evolution of an open system caused by its coupling to vacuum may 
become potentially observable. Here, we show that the modification of vacuum fluctuations 
induced by the presence of boundaries provides such a possibility to unveiling quantum vac- 
uum fluctuations via geometric phase. At this point, it is worth noting that geometric phase 
has recently been proposed as a possible way to detect the Unruh effect in Ref. [25] and 
later in Ref. [26]. 

The system we study contains an inertial two-level atom in interaction with a bath of 
fluctuating quantum electromagnetic fields in vacuum at a fixed distance to a reflecting 
boundary. The Hamiltonian of the whole system takes the form H = + + H' . Here 
Hg is the Hamiltonian of the atom, which, for simplicity, is taken to be = ^hLOQa^, 
where is the Pauli matrix, and uq is the energy-level spacing of the atom. Hrp is the 
Hamiltonian of the electromagnetic field, of which the explicit form is not relevant here. In 
the multipolar coupling scheme [27], the interaction Hamiltonian H' takes the form H'{t) = 
— er ■ E(x(r)) = —^Ylmn^mn ■ E(x(r))(Tmn) where e is the electron electric charge, er the 
atomic electric dipole moment, and E(a;) the electric-field strength. Here, the dipole moment 
must be kept fixed with respect to the proper frame of reference of the atom; otherwise, the 
rotation of the dipole moment will bring in extra time dependence in addition to the intrinsic 
time evolution [28]. Since neither r nor E(x) is a world vector, the interaction Hamiltonian 
H' is ambiguous when we deal with the situation of moving atoms. However, we can write 
the interaction Hamiltonian H' in a coordinate invariant form as H' = —er^F^^{x)u'^, 
where F^j, is the field strength, is a four-vector such that its temporal component in 
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the frame of the atom vanishes and its spatial components in the same frame are given by 
r, and u'^ is the four velocity of the atom. Since we choose to work in the frame of the 
atom, u" = (c, 0,0,0), and this coordinate invariant interaction Hamiltonian reduces to the 
form given above [28, 29]. The dynamical evolution of the two-level-atom subsystem will be 
studied in the paradigm of open quantum systems. Here, let us note that the theory of open 
quantum systems has been fruitfully applied to understand, from a perspective different 
from the traditional, the Unruh, Hawking and Gibbons- Hawking effects, in Refs. [30], [31], 
and [32], respectively. 

The initial state of the whole system is characterized by the total density matrix ptot = 
p(0) (S) |0)(0|, in which p(0) is the initial reduced density matrix of the atom, and |0) is the 
vacuum state of the field. In the frame of the atom, the evolution in the proper time r of 
the total density matrix ptot satisfies 

dptotir) _ i 

— q;^ — - -p^,Ptot{r)\ . [I) 

We assume that the interaction between the atom and the field is weak. In the limit of 

weak coupling, the evolution of the reduced density matrix p(r) can be written in the 

Kossakowski-Lindblad form [30, 33-35] 

c?p(r) i 
dT ^ ~h 

where 



[i/eff, P(r)] + C{p{r)] , (2) 



^[p] = ^ ^ aij [2 ajp (Ti - (Ti(Tj p- p (Ti(Tj] . (3) 

The matrix ajj and the effective Hamiltonian i^cfr are determined by the Fourier and Hilbert 
transforms of the correlation functions, 
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G^x -y) = ^Yl {Omx)E,{ym , (4) 



which are defined as follows: 



/oo 
rfre*^^G+(x(r)) , (5) 
'OO 

iri .1^ uj - X 



Then the coefficients of the Kossakowski matrix aij can be expressed as 

ttij = A5ij - iBeijkSks - ^^is^js , (7) 



in which 



A = -^[g{uj,) + g{-ujo)] , B = -^[g{ujo) - q{-ujo)] . (8) 



The effective Hamiltonian ifefr contains a correction term, the so-called Lamb shift, and one 
can show that it is given by replacing uq in Hg with a renormalized energy-level spacing Vt 
as follows [30]: 

1 hi 
iJcff = 2 ^^^3 = 2(^0 + -[/C(-a;o) - K.{ujq)]] . (9) 

Assuming the initial state of the atom is |^/'(0)) = cos ||+) + sin ||— ), one can show that 
the time-dependent reduced density matrix of the atom is given by 

p-AAr 2 9 I B~A ( -AAt _ i \ 1 -2Ar-inT • n 

P{t)= I 2+ 2A 2^ Sin^ I ^^Q) 

ie-2^-+*^-sin^ l-e-4^-cos2|-^(e-4^--l) 

which evolves nonunitarily. The geometric phase for a mixed state under a nonunitary 
evolution can be defined as [15] 

7, = arg l^fj VA.(O)A,.(T)(0,(O)|</.,(T))e-/o"(<A'=MI<^'=M)'^-j , (11) 

where Afc(r) and |0a:(t)) are the eigenvalues and eigenvectors of the reduced density matrix 
p(r). In order to find the geometric phase, we first calculate the eigenvalues of the density 
matrix (10) to get A±(r) = |(1±?7) , in which 77 = \/ p1 + e~'^^'^ sin^ 6 and ps = e"^"^"^ cos 6* -f- 
^(e"^"^"^ ~ !)• It is easy to see that A_(0) = 0. As a result, contribution only comes from 
the eigenvector corresponding to A+ , 

l0+(r))=sin||+) + cos|e^^n-), (12) 

where 



tan^ = ,/^. (13) 
2 V^-P3 ^ ' 



The geometric phase can then be calculated directly using Eq. (11), 



7, = -fi /"^cos^^dr. (14) 







Now, we calculate the geometric phase of an atom in the vicinity of a reflecting boundary. 
To do so, we need the two point functions for the electric fields, which can be found from 
those of the four-potentials, 

D^'^(x,x') = {{)\A^'{x)A''{x')\{])=D^^\x,x') + D''^\x,x') , (15) 



in which Dq^{x,x') is the two point function in the Minkowski vacuum without boundaries, 
and D^'^{x, x') is the correction induced by the presence of the boundary which can be 
calculated using the method of images. In the Feynman gauge, we have, at a distance z 
from the boundary, 

nf^^fj. = _1 '1 (TR) 

^ 47r%oC [{ct - ct' - leY - (x - x'Y - {y - y'f - [z - z' f] ' ^ ' 

h y-^'^ ) 4^2^gc [{ct - ct' - tey - (x - x'y -{y- y'Y - {z + z'f] ' ^ ' 

where rj^^^ = diag{l, —1, —1, —1), = (0, 0, 0, 1), and e — ?■ +0. The electric field two-point 
functions can be expressed as a sum of the Minkowski vacuum term and a correction term 
due to the boundary: 

{EMr))EM'^'))) = {EMr))EMr')))o + {E,{x{T))E,{x{r'))), , (18) 

where 



(0|E,(x(r))i?,(x(r'))|0)o = ^^(W., - 



he 

^ (x - x'y + {y- y'Y + (z - z'Y - {ct - at' - ley ' ^^^^ 



he 

{0\Ei{x{r))Ej{x{r')mb = --^m,-2n,nj)dod',-d,d'^] 

Ait Sq 

1 . s 

X f20l 

{x - x'f + {y- y'Y + {z + z'Y - {ct - ct' - ief ' ^ ' 

Here d' denotes the differentiation with respect to x' . 

Let us now consider an atom moving in the x-direction with a constant velocity v at a 
distance z from the plane, so the trajectory is given by 

t(r) = 7r , x(r) = Xq + inr , ?/(r) = y^ , z{t) = z , (21) 

where 7 = (1— f ^/c^)~^. Here let us note that the electric-field two-point functions Eqs. (18)- 
(20), and the trajectory Eq. (21) are described in the laboratory frame. Since the evolution 
of the atom is studied in the frame of the atom, a Lorentz transformation is required to 
get the electric field two-point functions in the proper frame of the atom from Eq. (18) to 
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Eq. (20): 

(0|E.(x(r))E,(x(r'))|0)o 



he 



6. 



vr^^o (c Ar — ie)^ 
|E.(x(r))E,.(x(r'))|0), = (0|E,(x(r))E,(x(r'))|0), 

|E,(x(r))E,(x(r'))|0), ^ 



TT^So [ (c Ar — iey — Az'^f 



vr^eo [ (c Ar - iey - Az'^f ' 
The correlation function and its Fourier transform can be calculated as follows: 

o2 



(22) 



where 



fx{>^,z) = fy{\z) 



e^\{-\n\+)?X ' 



3c3 



" 8A323 



2Xz 2Xz AX^z^ , 2Xz 
cos h ( — ^ 1 ) sin 



c 



c 



2Xz 2Xz 2Xz 
cos sin 



(23) 
(24) 

(25) 
(26) 



c c c 

and ^(A) is the standard step function. Thus the coefficients of the Kossakowski matrix a^j 
and the effective level spacing of the atom can be written as 



n = uo + -K^TP r duu'a,{l - f.iuo, z)) (— ^—] 



(27) 



(28) 



where 70 = e^|(— |r| 



l/Sireohc^ is the spontaneous emission rate in vacuum without 



boundaries, and = |(— |r|+)P • Then the geometric phase can be found using 
Eq. (14), 
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ill 



1 - e^^^ + cos 9 



Qdr . 



(29) 



V^e^^^ sin^ ^ + (1 - e^^^ + cos Oy , 
So, the phase accumulates as the system evolves, although the accumulation with time is 
not linear as in the unitary evolution case. For a single period of evolution, the result of 
this integral can be analytically expressed as 

Q 



7, = -(F(27r)-F(0)) 



(30) 
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where function F{(p) is defined as 




sgn(g) In (2g2e-4AW-o -Q^ + 2Q S{ip) 6-^^^/"°)) , 



(31) 



in which S{(f) = ^ e^A^I^<^ - (AA^I^i,Qi + Q2 , Q = i + cos6', and sgn((5) is the standard 
sign function. 

In order to examine the behaviors of this phase, we perform, for small 7o/wo; which is 
true in our current discussions as we will see later, a series expansion of the geometric phase 
for a single quasi-cycle and find, to the first order, ^ 



The first term — 7r(l — cos &) in the above equation is what we would have obtained if the sys- 
tem were isolated from the environment, i.e., a bath of fluctuating vacuum electromagnetic 
fields and the second term is the correction induced by the interaction between the atom 
and the environment. Here fi{u}o, z) are oscillating functions of distance z with a position- 
dependent amplitude. For an atom polarized in an arbitrary direction, the polarizations of 
the atom in the tangential directions and in the normal direction of the boundary contribute 
differently to the correction of the geometric phase. If the atom is polarized in the tangential 
direction, as the atom approaches the boundary [z — )■ 0), the correction of the geometric 
phase vanishes, since fx{ojQ,z) and fy{uo,z) approach zero, which can be attributed to the 
fact that the tangential components of the electric field vanish on the conducting plane. 
However, if the atom is polarized in the normal direction, fz{u}o, z) —1 as z 0, and the 
correction of the geometric phase is twice that of the free space case. This can be under- 
stood as the fact that the reflection at the boundary doubles the normal component of the 
fluctuating electric field. When the distance z approaches infinity, the modulation functions 
fi{u}o, z) approach zero, and the result reduces to that of the unbounded Minkowski vacuum 
case. So, due to the modification of the vacuum fluctuations caused by the reflecting plane, 
the vacuum fluctuation induced geometric phase becomes position dependent. Now let us 
estimate how large the phase difference is. If we assume that |(— |r|+)| is of the order of 

^ Here we have omitted the Lamb shift terms, since it is obvious that these terms eontain a faetor 70/wo 
and they wiU only contribute to the phase at the second and higher orders of 70 /wq ■ 



7r(l — cos 9) — TV 



2 7o 
2^0 



^ai(l-/i(a;o,^))(2 + cos^) sin^ ^ . 



(32) 
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the Bohr radius oq, and uq is of the order of Eq/H, where £"0 = —e^ /Sneocio is the energy of 
the ground-state, then '-fo/uo is of the order of 10~^. For a fixed z, the environment induced 
geometric phase (the second part of Eq. (32)) reaches its maximum when 9 ^ 1.354, which 
is in the vicinity of ^ = it/2, i.e., an equal superposition between the ground and excited 
state. (The numerical results below are based on 6* = 7r/2.) 

Based upon the discussions above, an experiment that aims a direct detection of vacuum 
fluctuations can be in principle designed. One first prepares two-level atoms in a superposi- 
tion of upper and lower states in a Ramsey zone. The atoms are then set in two paths which 
are both parallel to a reflecting plane but with different distances from it. After a certain 
time of evolution, we let atoms from different paths meet by certain means, for example, by 
applying a laser pulse to atoms from one path to change its direction of motion, and take an 
interferometric measurement. For a practical experimental implementation, we want, on one 
hand, the difference in the distances of the two paths to the boundary to be large enough so 
as to generate an appreciable phase variance. On the other hand, however, we also want this 
difference to be negligibly small as compared to the length of the paths parallel to the plane, 
so that the phase accumulated due to the vertical motion of the atoms before they meet 
can be neglected. A compromise can be achieved for atoms whose transition frequencies are 
in the microwave regime, for example, Uq ~ 10^ s^^, which is physically accessible [36]. If 
the distances to the plane of the parallel paths are chosen respectively as ~ 10~^ m and 
~ 10~^ m, then one can show by integrating Eq. (29) that the geometric phase variation can 
reach ~ 10^^ rad for an evolution time of ~ 10"'^ s. In current cold atom interferometric 
experiments, the speed of the atoms is ~ 1 m/s, so the time the atom moves vertically 
is ~ 10~^ s, which is two orders of magnitude less than the parallel evolution time, and 
thus the phase accumulated during this period can be neglected. Here we emphasize that 
the geometric phase is independent of the velocity of the atom. This can be seen from the 
Fourier transform of the electric-field correlation functions Eqs. (24)-(26), which determine 
the coefficients of the dissipator A and B (Eq. (8)), and then the geometric phase (Eq. (29)). 
We specify the velocity of the atom here only to ensure the geometric phase generated by 
the motion in the parallel direction will dominate. In Fig. (1), we plot the phase difference 
between one atom whose trajectory is fixed at zq = 1 fim, and the other that varies from 
z = 1 fim to z = 100 yum, which shows that the phase difference increases with the distance 
between the two atoms monotonously. Now, we estimate how the geometric phase would 
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change when the trajectories fluctuate by an amount bz. We assume one trajectory is fixed 
at 2o = 1 /^ni; and the other fiuctuates from z = 10 /im to 2; + 5z = 10.1 /im. The phase 
difference between the two cases is ~ 10"^ rad, which is two orders of magnitude smaller 
than 10"'^ rad. So the geometric phase is robust against small fluctuations of the distances 
5^, as long as biz is small compared with zq. Another effect that should be taken account 
of is that, in reality, a metal plate does not reflect electromagnetic waves completely. As a 
result, an excited atom also decays nonradiatively, i.e., the energy is not only transferred to 
the free space as photons but also to the absorbing metal as heat. The nonradiative decay 
rate takes the well-known form 7„on/7o = /? -2"'^ (see, e.g., Ref. [-37] and references therein, 
and Ref. [38] based on a fully canonical quantum theory). This would have an effect on the 
environment induced geometric phase as can be seen from Eq. (32), that is, to the flrst order, 
the correction is proportional to the spontaneous emission rate. For conductors, /3 is typi- 
cally of the order of ~ 10"^^ cm^ [37], and we are considering 2; ~ 1 /im, so 7non/7o ~ 10~^. 
Therefore, the contribution of the nonradiative decay to the phase can be neglected. We 
must point out, however, that a subtle issue actually exists in any practical implementation 
of our proposal, that is, a cancellation of dynamical phases that the atoms may acquire dur- 
ing the evolution, which is very tricky for systems under non-unitary evolutions like what we 
are considering here [39]. A possible alternative might be to determine the geometric phase 
directly in a tomographic manner by measuring elements of the reduced density matrix of 
the atom rather than to perform an interferometric experiment as what is actually pursued 
in [24] . 
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<5(rad) 




FIG. 1: Phase difference S as a function of z. Here z and zq are the distances of the trajectories of 
the two atoms to the plate respectively, with parameters zq = I /xm, wq = 3 x 10^ Hz, 7o/c<Jo = 10~^, 
and e = tt/2. 
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